Inhomogeneous superconducting state in mesoscopic thin-walled cylinders 
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We theoretically investigate the appearance of spatially modulated superconducting states in 
mesoscopic superconducting thin-wall cylinders in a magnetic field at low temperatures. Quanti- 
zation of the electron motion around the circumference of the cylinder leads to a discontinuous 
evolution of the spatial modulation of the superconducting order parameter along the transition 
line T C (H). We show that this discontinuity leads to the non-monotonic behavior of the specific 
heat jump at the onset of superconductivity as a function of temperature and field. We argue that 
this geometry provides an excellent opportunity to directly and unambiguosly detect distinctive 
signatures of the Fulde-Ferrell-Larkin-Ovchinnikov modulation of the superconducting order. 
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Mesoscopic systems both serve as a platform to inves- 
tigate fundamental quantum physics of solids and are a 
testing ground for potentially transformative future de- 
vices [iHSl- Of special interest in this context are in- 
teracting systems exhibiting interplay of the collective 
emergent properties with the quantum physics of single 
particles. Small superconducting samples of nontrivial 
topology provide an example of such interplay since the 
global phase of the pair condensate and the phases of 
single particle wave functions respond differently to the 
applied magnetic field. 

In mesoscopic superconducting rings of radius R ~ £oj 
where £o = vf/2ttT c (H = 0) is the superconducting co- 
herence length, and vf is the Fermi velocity, this leads to 
a doubling of the period of the oscillations of the transi- 
tion temperature, T c , as a function of the magnetic flux, 
<f>, through the ring, relative to the well-known Little- 
Parks (LP) effect [5|-tZ|. The small ring radius, R, im- 
plies that each single electron state can be labeled by its 
angular momentum, n, in units of h, and each particle 
acquires an additional phase due to the magnetic flux, 
$ = ttR 2 H when circling around the ring. In the ab- 
sence of the field the wave function of the Cooper pair 
has a net zero angular momentum, as the time-reversed 
states with ri\ = —ri2 form a bound state. In contrast, 
under the applied field the non-zero quantum number 
I = n\ + ri2 partially compensates the net flux and max- 
imizes the transition temperature, T c . Therefore, for a 
small ring 0], T c is a periodic function of <j) = <J>/(2<J> ) 
with the flux quantum $o = hc/2e, while for a large ring 
the periodicity is solely due to the flux experienced by an 
electron pair with charge 2e, i.e. T c is a periodic function 
of $/$ (LP effect). 

Under these assumptions, there is no overall suppres- 
sion of the maximal T c as the magnetic field increases: 
at integer values of <j> the orbital coupling of the field to 
the individual electrons can be fully compensated by the 
finite angular momentum I of the Cooper pair. The de- 
struction of superconductivity in such a geometry must 



then occur via paramagnetic (Zeeman) coupling of the 
electron spins to the field, which raises the energy of the 
singlet bound state: inclusion of this coupling is essen- 
tial for developing a complete picture. Consequently in 
this Letter we consider the combined effect of the or- 
bital and Zeeman effect on superconductivity, and an- 
alyze a mesoscopic thin-walled cylinder with the field 
along the axis. The cylinder geometry allows the forma- 
tion of a spatially- modulated [g, Q (Fulde-Ferrell-Larkin- 
Ovchinnikov, FFLO) superconducting state that enables 
pairing under high Zeeman field. We show below that 
a) this state occurs even if the cylinder is made out of 
materials where superconductivity is not paramagncti- 
cally limited in the bulk; b) signatures of such a state are 
much more prominent and easily identified in this geome- 
try with R ~ £o than either in bulk materials or fiat thin 
films, and therefore mesoscopic systems offer a unique 
chance to detect the FFLO state that has remained elu- 
sive for nearly 50 years since it was first predicted. 

Our main results are shown in Fig. [TJ While at low 
fields the variation of the transition temperature T C (H) 
is well described by the Little-Parks (LP) periodicity, 
at higher fields there is an extended region in which 
the superconducting order parameter is modulated along 
the cylinder axis. Near the phase boundary T C (H) in 
this regime, the superconducting phase exhibits alter- 
nating regions of phase-modulated (FF) and amplitude- 
modulated (LO) order; however, the LO phase becomes 
stable at lower T. The wavevector, q z , of this modu- 
lation has non-analytic dependence on H, due to the 
interplay between the finite size effects and the LP os- 
cillations. The heat capacity jump at T C (H) varies dra- 
matically along this sequence of transitions, enabling a 
direct identification of the modulated states. 

Although the possibility of FFLO states has been dis- 
cussed in bulk materials such as the heavy-fermion su- 
perconductor CeCoIns U 12 1 and organic superconduc- 
tors 13-1(3], it is difficult to design a "smoking gun" 



experiment that unequivocally points towards such a 
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FIG. 1: (Color online) Superconducting thin- wall cylinder 
with a small radius R in a magnetic field parallel to the cylin- 
der axis (z-axis). Sketch shows the spatial modulation of the 
order parameter in the LO phase. Right-panel: phase dia- 
gram in the H-T plane. Shaded regions correspond to the 
phase-modulated (FF) states. 

state. In real bulk systems both orbital and paramag- 
netic coupling suppress superconductivity, and inhomo- 
geneity arises due to both. The former effect leads to 
proliferation of the vortices. Each vortex contains ex- 
actly one flux quantum for the Cooper pairs, $0j which 
corresponds to a 27r-phase winding of the superconduct- 
ing order parameter around the vortex core. In that sense 
a thin-walled ring or cylinder can be viewed as a corc- 
less "supervortex" of multiple flux quanta, with a phase 
winding 2irl. In contrast, paramagnetic pairbreaking fa- 
vors FFLO states. Recall that in a singlet superconduc- 
tor in the absence of Zeeman splitting the Cooper pair 
comprises electrons in time-reversed states, which have 
equal energies, and therefore are unstable towards for- 
mation of a bound state. With paramagnetic coupling 
the states with opposite spins have equal energies if they 
have a field-dependent momentum mismatch g, and the 
modulation of the FFLO state originates from this finite 
center-of-mass momentum of the Cooper pairs. On a ring 
this is directly equivalent to the angular momentum mis- 
match, and hence there is no clear distinction between 
the vortex and the modulated phases. 

We consider a long hollow superconducting cylinder of 
radius R, which, in the absence of a magnetic field, is 
described by the Hamiltonian 

H = ]T£(p)c PiCT c p , CT -A^Bt (q) £ (q)5 

a,p q 

#(q) = \ (- icr yk/3 a -p+a, Q c p+ a i/3 ■ (1) 
p,a,/3 

Here c Pja is the annihilation operator for an electron with 
momentum p and spin projection a, A is the strength 
of the pairing interaction, and we assumed singlet s- 
wave superconductivity. For small R the motion around 
the circumference of the cylinder is quantized, while the 
momentum along the axis (z) is continuous, so that 
p = {m/R,p z ) with m integer, and the quasiparticle en- 
ergy takes the form 



where M is the electron mass and fi is a chemical po- 
tential. Here the summation over the momenta means 

E P =[(2 7 r) 2 i?]" 1 E me z/^- 

The magnetic field threading the cylinder along its axis 

leads to a Zeeman splitting of the single particle energy 

levels by H z = J2a. P ° h 4,^ S p^' wnere h = Mb .g# , fi B 
is the Bohr magneton, and g is the g-factor of the quasi- 
particles in the crystal. At the same time the momen- 
tum operator has to be replaced by its gauge-invariant 
counterpart, p — > p + \e\A. For our model of a thin- 
walled cylinder the vector potential A = — (pH R/2, and 
(p is the azimuthal angle around the cylinder. Hence 
|A| = const on the cylinder. 

In the superconducting (SC) phase the pair field 
A(q) = X(B(q)), where (...) denotes the thermal av- 
erage, acquires a non-zero value. Due to the cylindrical 
geometry, A(r) can be expanded in the Fourier series, 

A(r) H A | e^e 4 *'*. ( 3 ) 

The uniform SC state at H = only has Co,o 7^ 0, while 
a single component C^o with flux-dependent I ^ char- 
acterizes the SC transition under orbital coupling to the 
field and gives rise to the Little-Parks effect. If the cylin- 
der were to unfold into a two-dimensional (2D) plane, I 
and q z would become components of a 2D vector q, and 
in response to a Zeeman field a state with q ^ would 
be realized (FFLO state). The cylindrical geometry is 
unique since / 7^ gives the magnetic flux through the 
cylinder, and therefore the additional FFLO modulation 
is only along the axis, q z 7^ 0. Thus, this setting has 
an advantage over other ways to achieve FFLO states 
(such as bulk paramagnetically-limited superconductors, 
imbalanced fermionic atomic gases or thin SC films), in 
which the direction of the modulation is arbitrary, and 
therefore complex states with simultaneous modulations 
along different directions may be favored (T7I-I22I] . 

To study the transition between the normal and 
SC states we use the Ginzburg-Landau (GL) expan- 
sion of the free energy, J" G l = (I, q z , T, H )| A | 2 + 
a^(l,q z ,T,H)\A \ 4 , At each H,l,q z , the temperature T 
where aA 2 ) (/, q z , T, H) becomes negative (if it exists) indi- 
cates a putative second order transition from the normal 
into the SC state with given values of l,q z . The high- 
est of these temperatures is the physical transition point 
T C (H) into a state with the corresponding I, q z . The nec- 
essary condition for the continuous transition is that the 
quartic coefficient a^ 4 \l, q z , T, H) remains positive at the 
transition point. 

We evaluate the coefficients a/ 2 ' and a' 4 ) using 
the Green's function formalism for the Hamiltonian, 
Eq. |T]) with the Zeeman and orbital coupling terms. 
We use a quasi-classical approximation for the nor- 
mal state Green's function — (T T ip a (r, t) ip^(r', 0)) ~ 

TJ2 £n z~ EnT Qe n A v ~ r ') eM ds ' A(s) , where the inte- 
gral in the exponent is evaluated along a straight line. 
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This approximation smears out the even-odd flux peri- 
odicity for a ID ring, but this periodicity is already bro- 
ken by Zeeman coupling, and hence the approximation is 
adequate for our goals. We obtain for the quadratic term 

a( 2 )|A | 2 = J drA*(r)(~-|££#( e „,<7))A(r), 

P 

K(e n ,a) = g ent<7 (p)g. £n ,. a (-p + Il), (4) 
where Q ena (p) = (ie n — £ l (p)+(7h)' 1 is the Fourier trans- 
I 



form of C? eiijCr (r), e n is a fermionic Matsubara frequency, 
and II = -i V + 2|e|A. 

The effects due to small ring size i? ^ £o are con- 
tained in the discrete sum over integers ra in V . Wc 

rjp ^P 

use the Poisson summation formula [7[, X^mez 3( x ~ m ) = 
J2kez e l2vkx to elucidate these effects: fc = gives the 
continuum result for a 2D superconductor, and higher or- 
der terms, k > 1, account for the finite size corrections. 
After a straightforward calculation, we find 



,(2: 



M 
2^ 



In 



2 [ ~9^ COS { 2lT k RpF cos (^p)) H 



fe>0 



n>0 



^ Jo 2tt 

™+5 



(5) 



4 ^ n 

S £ . (7 



e~S'fc 7r(2 ™ +1)fc|cos(¥ ' p)l 



where we defined the product 
v F Q = 2irT c (£ q z sm(ip p ) + 



$ 



R 



■cos(p p )). (6) 



Here, ip(z) is the digamma function and the Fermi mo- 
mentum is pf = y/2MfI. In the second line of Eq. (|5|) we 
neglected terms of order T 2 /^ 2 . Because of the exponen- 
tial decay of the terms with increasing fc > in the last 
term of Eq. ([5]), below we keep only the first finite size 
correction, fc = 1. Wc checked that incorporating higher 
fc does not qualitatively change our results. 

Fig.[5]shows the upper critical field and the parameters 
l,q z of the modulation of the superconducting order pa- 
rameter at transition. Hereafter we consider £,oPf = 100, 
and R ~ 3£o- In the cylindrical geometry the inhomoge- 
neous superconducting states emerge even for the mate- 
rials that do not support FFLO modulation in the bulk: 
we present the results for the paramagnetic parameter 
a M = .9/ i s c E>o/(7rCo ^c) = 0.6, which corresponds to the 
Pauli limiting field Hp « 4i7°| b , so that the bulk ma- 
terial is a conventional orbital-limited typc-II supercon- 
ductor. The scalloped shape of the boundary of the su- 
perconducting region is the consequence of the LP effect, 
and the overall suppression of T c with increased field is 
due to the paramagnetic pairbreaking. Below a charac- 
teristic temperature, which is non- universal and different 
from the T* = 0.56T c o for bulk Pauli limited supercon- 
ductors, the inhomogeneous pairing along the cylinder 
axis (q z ^ 0) becomes advantageous, and the FFLO state 
appears. 

Fig.[2]shows that the modulation wave vector, q z , along 
the transition line exhibits a "sawtooth" pattern, quite 
distinct from the uniform increase in q in the standard 
picture of Pauli-limitcd superconductors. This feature is 



due to the effective discretization of the modulation in 
Eq. ^ . For a 2D sheet the role of the winding number 
I is taken by a continuous variable q x , and it is the net 
9 = \/<7 2 + <7 2 that ensures matching of the energy of the 
two electrons in a Cooper pair with the center of mass 
momentum q. In contrast, in the cylindrical geometry the 
choice of I is determined by the net flux, $, and therefore 
the momentum q z adjusts to this selection, and exhibits 
discontinuities at the points where transitions between 
winding numbers I and / + 1 occur. The detailed balance 
between q z and I depends on the finite size quantum cor- 
rection term, second line in Eq. (|5|) . Note that the prefac- 
tor of the first non- vanishing fc term, cos[2irRpp cos(y p )], 
has the same angle-dependence in the momentum space 
as the the LP term, \l — $/<&o| cos(ip p )/R, and is out 
of phase with the FFLO modulation that enters with 
sm((p p ) in Eq. (|6]). Consequently the details of the switch- 
ing between values of I and q z depend on the value of RpF 
(close to integer vs. half- integer), but the qualitative pic- 
ture remains unchanged. 

The discontinuous behavior of the modulation q z (T) 
is reflected in the experimental properties that allow un- 
ambiguous determination of the modulated state. Fig. [3] 
shows the specific heat jump at the superconducting tran- 
sition for different fields. We verified that the quartic 
term, cS 4 ^ (T C (H) , I , q z ) , remains positive along the en- 
tire transition line, and therefore the transition is al- 
ways of second order. In that case whether the phase- 
modulated (FF) or the amplitude-modulated (LO) state 
is favored is simply determined by comparing the mag- 
nitude of aW(T c (iJ)) for the FF and the LO states: 
smaller value corresponds to the greater condensation 
energy and a more stable phase. We find that in the 
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vicinity of the discontinuous drop of the modulation q z 
the FF state is favored, and is superseded by the LO 
state as q z increases within the realm of each fixed wind- 
ing number I. The specific heat jump at the transition 
is given by AC/T C (H) = ([a®]') 2 /2a< 4 > (T C (H)), where 
[aW]' = {da ( - 2 \T,H)/dT) T=Tc(H) . The heat capacity 
exhibits significant enhancement on transitions between 
different winding numbers. The enhancement of the spe- 
cific heat jump is measurable (and should be much more 
prominent in the hollow cylinder geometry, where it can 
be detected, for example, by the ac calorimetry tech- 
nique), and therefore can serve as experimental proof of 
the existence of the FFLO-like modulations of the super- 
conducting order in mesoscopic cylinders. 

It is well known that the modulated state is sensitive to 
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FIG. 2: (Color online) Structure of the modulated state for 
R/£o = 3 and a = 0.6. (a) The upper critical fields H C 2 nor- 
malized by <3>o/7r£o- Inset and dashed lines show H C 2(T) for 
a given angular momentum I pairing state as indicated. Solid 
symbols denote the physical transition. Note the switching 
between the FF (triangles) and LO (circles) states along the 
transition line, (b) FFLO modulation wave vector q z nor- 
malized by l/£o for each I (dotted lines), and at the physical 
transition as in panel (a). The circles (triangles) denote the 
stability regions of the LO (FF) state. Note the non-analytic 
behavior of q z exhibiting kinks and discontinuous jumps at 
temperatures denotes by arrows. 
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FIG. 3: (Color online) The specific heat jump at the onset of 
the SC order, AC (T C (H))/T C (H). The notations are the same 
as in Fig(2] The inset shows the evolution of the specific heat 
near the discontinuity in the modulation q~ at T/T c = 0.29. 
All quantities are normalized by their values at T C (H — 0). 



impurities [231 ] . and it is likely that in the experimental 



realization of the proposed geometry the superconduc- 
tor will be moderately disordered. We checked that ac- 
counting for impurity scattering, as long as the transport 
remains ballistic and not diffusive, does not lead to the 
complete destruction of the modulated states, and pre- 
serves the non-monotonic behavior in the q z (T) curve. 
The situation is similar to that found in the detailed 
treatment using the T-matrix approach, where the FFLO 
bulk state remained stable over a substantial range of the 
impurity concentrations 24, 25j |. 

To conclude, we find novel spatially-inhomogcneous 
superconducting states, exhibiting both the Little-Parks 
and FFLO phenomenology, can emerge due to the vector 
potential and Zccman coupling induced by a magnetic 
field threading a thin cylinder. Our principal motiva- 
tion was conventional superconductors, for which the co- 
herence length can be well in excess of fOOOA. In this 
setting, the relevant sample sizes are experimentally ac- 
cessible and we believe that the predicted variations in 
the specific heat jump can be found under realistic condi- 
tions, providing a possible "smoking-gun" experiment for 
detecting the FFLO state. In principle, recent achieve- 
ment of artificial electroma gnet ic vector potentials in 
cold atomic gas experiments [261 ] . cold gases confined to 



a cylindrical trapping geometry could provide another 
possible setting for realization of similar phase diagrams 
(with an imposed population imbalance playing the role 
of the Zeeman coupling [27|). However, given the dif- 
ficulty of directly measuring the heat-capacity jump in 
a trapped cold atomic gas, it may currently be easier to 
detect the effects in small-size superconducting systems. 
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